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Abstract 

We show here for the one-dimensional spin-1/2 ANNNI (axial-next-to-nearest- 
neighbor-Ising) model in an external magnetic field that the linear density of Yang- 
Lee zeros may diverge with critical exponent a = —2/3 at the Yang-Lee edge singu- 
larity. The necessary condition for this unusual behavior is the triple degeneracy of 
the transfer matrix eigenvalues. If this condition is absent we have the usual value 
a = —1/2. Analogous results have been found in the literature in the spin-1 Blume- 
Emery-Griffths model and in the three-state Potts model in a magnetic field with 
two complex components. Our results support the universality of o" = — 2/3 which 
might be a one-dimensional footprint of a tricritical version of the Yang-Lee-Edge 
singularity possibly present also in higher-dimensional spin models. 
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1 Introduction 



In spin-1/2 models in an external magnetic field if, the partition function for finite num- 
ber of spins is proportional to a polynomial on the variable u = exp(— 2if/fcT). Since 
polynomials are basically specified by their zeros, those zeros furnish all relevant physical 
information. Given that the partition function is a sum of exponentials (positive num- 
bers) with positive coefficients, it is clear that we can only have zeros in the u- variable for 
complex magnetic fields. Such zeros for complex magnetic fields are called Yang-Lee zeros 
(YLZ) and are naturally studied on the complex-u plane as we do here. Their relevance 
in the study of phase transitions has been pointed out in 1952 by C. N. Yang and T. D. 
Lee, see [1]. 

In several spin models the YLZ tend, in the thermodynamic limit, to form continuous 
curves on the complex u-plane. In particular, it has been rigourously proved in [2J that 
all YLZ of the spin-1/2 Ising model, even before the thermodynamic limit, lie on the unit 
circle |m| = 1. This circle theorem has been generalized, for example, to higher-spins Ising 
models |3] and to include other interacting terms |4j. As a rule, the zeros tend to leave 
the unit circle as the ferromagnetic (FM) Ising coupling J becomes smaller as compared 
to other couplings. For more references see the review work [5]. At T <Tc the YLZ tend 
to pinch the positive real axis on the complex-u plane at the first-order and second-order 
phase transition points as we approach the thermodynamic limit. However, ii T > Tc 
they accumulate at the endpoints of the edges {ue = Gxp{—2f3HE)) of the curves with a 
divergent linear density p{u) ~ |m — ueI"^ with a < 0. 

As first noticed in [6] for the two-dimensional Ising model, the power-like behavior is 
independent of the temperature as long as T > T^. The universality of the exponent a 
was explained in |7] as a result of a usual critical point described by a field theory with 
i interaction vertex. The corresponding endpoints ue have been called Yang-Lee edge 
singularities (YLES). In D = 2 dimensions, the use of conformal field theory |H] predicts 
a = —1/6. This result has been verified for the Ising Model numerically [9l [10] and also 
experimentally from magnetization data pT] . 

According to arguments given in [12] one should have a = —1/2 in one dimension. 
Notwithstanding, even in D = 1 the exact position and density of the YLZ is not known 
analytically. One exception is the one-dimensional spin-1/2 Ising model. In this case the 
linear density was already known exactly [2] furnishing a = —1/2. Numerical works have 
confirmed a = —1/2 in several one-dimensional spin models fl3 \ [T ^ [T5 | [T6] including, see 
[H], the same model discussed here. An early exception is a special type of three-state 
Potts model [TH]. In such model the spins are coupled with a magnetic field with two 
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complex components. By fine-tuning the couplings of the model it has been obtained 
another type of YLES with o = —2/3. More recently, inspired by the work of we 
[T9| [20] have shown that a = —2/3 is not a peculiar feature of the special Potts model 
used in [18] but it is also present in the more familiar spin-1 Blume-Capel |2T] and spin-1 
Blume- Emery- Griff ths (BEG) [22j models. Once again a fine-tuning of the couplings is 
required, which is typical of a tricritical phenomenon. In particular, we need to have a 
triple degeneracy of the transfer matrix eigenvalues. So far one has found this unusual 
value for a only in spin models with three states per site. Here we investigate the one- 
dimensional spin 1/2 axial- next-to-nearest-neighbor Ising (ANNNI) model, and confirm 
that the same unusual critical behavior also appears in models with two states per site and 
next-to-nearest-neighbor interaction. Our setup is based on the transfer matrix solution 
and the use of finite size scaling relations which are shown to be satisfied by the YLZ 
close to the unusual YLES. 

2 The one-dimensional ANNNI model 

The ANNNI model was introduced in [23], see also p3], as a simple model to describe 
spatially modulated periodic structures observed in magnetic and ferroelectric materials. 
In its higher-dimensional {D > 1) versions the model has interesting physical applications, 
see [25] for a review work. We concentrate here in its one- dimensional version as a simpler 
laboratory to investigate the existence of other types of Yang-Lee-Edge singularities. The 
energy and the partition function of the spin 1/2 one-dimensional ANNNI model in an 
external magnetic field are given by: 

i i i 

Z = J2e-'' (2) 

where 5*, = ±1, J and K are coupling constants between nearest and next to nearest 
neighbor spins respectively and H is the magnetic field. Usually the model is defined 
with J > (ferromagnetic or FM) and K < (anti-ferromagnetic or AFM) couplings. 
Here we start from a more general standpoint where J and K are arbitrary real numbers 
and the magnetic field H can assume complex values. We use the notation: 

c = e-^-^, d = e^^, u = e^^^, A = m + 1/u = 2 cosh(2/3i/). (3) 
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The temperature is given in units of the nearest neighbor couphng, i.e., the range 
< T < oo corresponds toO<c<lifJ>Oorc>lifJ<0. Likewise — oo < K < 
imphes < d < 1 while < K < oo leads to d > 1. Here we use periodic boundary 
conditions, 5*^ = Si+N- The partition function can be found via transfer matrix [1] 



= Tr = Af + + + Af (4) 

where the Ai, A2, A3 and A4 are the eigenvalues of the transfer matrix of the model and 
they can be determined by the characteristic equation: 



A"^ - flsA^ + a2A^ - fliA + Oo = (5) 

where 



«o = -{j.-d'X (6) 



ai = dc\^—-d^] VA + 2 (7) 
'^-VaT2 (9) 



as 



c 



The Z2 symmetry {H — )■ —H) of Zjy is explicit in the factor y/A + 2 = ^ + e~^ ^ present 
in ai and a^. Notice that for d = 1 we recover the spin- 1/2 Ising model with only two 
eigenvalues. 

3 a = —1/2 and a = —2/3 (Analytic approach) 

The partition function Zn is proportional to a polynomial of degree in the "fugacity" 
u (lattice gas interpretation). All relevant information about Zn is contained in its zeros 
ZN{uk) = 0, = 1, 2, A^. Due to the Z2 symmetry one half of the zeros is the inverse 
of the other half. For large number of spins we assume, see [26l IT4j . that the zeros are 
determined by imposing that at least two of the eigenvalues of the transfer matrix have 
the same absolute value which must be larger than the other two remaining ones, i.e.. 
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(10) 



IA2I = |Ai| > |A,|, z = 3,4. (11) 

For large N the contributions of A3 and A4 can be disregarded and the partition function 
becomes 



Zj^^\^{l + e'^^) (12) 

Therefore the zeros are given by 



ypfc = (2A; - 1) -, k = l,2,..,N. (13) 

The equations ([5]) and ( ITOj) imply 

as = Ai(l + e^'^)+A3 + A4 (14) 

as = A^e*^ + Ai (A3 + A4) (1 + e*^) + A3A4 (15) 

ai = A? (A3 + A4) e*'^ + A1A3A4 (1 + e*'^) (16) 

ao = A?A3A4e^'^ (17) 



In ([I^-(|T71) we have four equations and four variables Ai, A3, A4, ip, ehminating Ai, A3 and 
A4 we find an equation which depends only on ip: 



F{A, c, d, (f) = 44 — 7 + 34 ao 0,2 — 28 ala2 — 2 afa^ + 6 ag ag — 40 ai a3 + 6 as a3 

— 24 ao ai a2 as + 34 ag as a^ + a^^ ag ag — 2 ao ag ag — 2 a^ ag + 6 ao ai as 

— 7a^a^ + 2 {40 a|^-a^ (24 a^ + 34 ai a3 - 29 as a^ + 6 a^) 

+ ao [al (29 as — 2 al) + al (4 as — a^) + ai as a^ (—19 as + 5 al)~\ 

— al [6 a^ + a2 + ai (—5 as a3 + ag)] } cos(v5) + 2 |31 al + al (—3 ai + as a3) 
+ ao [al + al (I8 as — al) + ai as ag (—8 as + ag)] — a^ (I6 a2 — 18 as ag 

+ 3 (8 ai a3 + ag)] } cos(2y9) + (40 ag — 2 a^ +14 ao al as — 16 a^ a^ — 26 a^ ai a^ 

— 2 ao ai al ag + 14 a^ as ag — 2 ag ag ) cos(3(y9) 

+ (20 aQ + 2 ao a^ as — 4 aQ a2 — 8 a^ ai a3 + 2 a^ as ag) cos(4y9) 

+ (8a|^-2aoaia3)cos(5v3) + 2ai^ cos(6v3) = (18) 
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A similai0 expression has been derived before in [T7] and analogous formulae for the one- 
dimensional spin-1 Blume-Capel and Blume-Emery-Griffiths models have appeared in [12] 
and |T9] respectiveljj^ . We interpret f|T8l) as a cubic equation for A = A{(p) such that when 
we plug it back in ([5]) we get two eigenvalues with the same absolute value according to 
f lTU]) . Equation flTS]) does not imply automatically the second condition f lTT]) . In practice 
we have to verify whether f lTT]) holds for each of the three solutions Ai{(p) , z = 1, 2, 3 of 
( !T8|) . At this point it is important to remark that, as its counterparts in [161 [13 ttS], 
the equation (fT8|) is symmetric under (f — )■ —(f. This symmetry is not accidental. It is 
a consequence of the permutation symmetry Ai A2 hidden in (fT4l) -( fT7l) . It becomes 
explicit when we perform {<f , Ai) — > {—ip , e^'^Ai). It is clear that after eliminating Ai , A3 
and A4 the resulting expression should be invariant under ip — )■ — Such symmetry will 
play a key role in determining the points in the parameters space of the model where the 
unusual critical behavior a = —2/3 shows up. Next we show how f lTSj) can be combined 
with a finite size scaling relation to determine the exponent a analytically. 

The closest zero Ui[N] to the YLES Ui[N — )■ 00], for large N, should obey [27] the 
finite size scaling relation: 



A.,(iV).„,|L]-„.[oo] = -^ = ^ . (19) 

Where i/h is the magnetic scaling exponent related to a via cr = (D — yh)/yh. The constant 
Ci is independent on the number of the spins N and = = L. 

It is known [281 EH] that the degeneracy of the two largest eigenvalues (Ai = A2) of the 
transfer matrix signalizes a second-order phase transition. Therefore, the YLES occurs 
at = 0. Thus, for large A^, we can assume that the closest zero Mi[A^] is obtained 
from the smallest angle as Mi[A^] = u[!f = (fi = tt/N], see ( |T3l) . We obtain ^(v^i) from 
the appropriate solution A{ipi) of ( |T8|) via u + 1/u = A. Expanding the result about 
A^ — 00, which amounts to know A{(f) in the vicinity of ip = 0, and comparing with ( ITQ]) 
we determine yh and a analytically. 

Although the exact solutions of the cubic equation ( ITSl) are cumbersome, they can be 
used to show that an expansion of A{ip) about p = contains only positive integer powers 
of p. The key ingredient is that the coefficients of the cubic equation (JT8l) are analytic 
functions of p. Therefore we can write down the large A^ expansion 

^The corresponding expression printed in [17], differently from ours, does not lead to a double degen- 
eracy of the transfer matrix eigenvalues at ip = as expected. 

^For the one-dimensional spin 1/2 Ising model the analogous of ([T8| is simply A — 2\J\ — cos {•p/'i). 
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A(0) + (^ 



dA 

dip 



+ 



d'A 



tfi=0 



, , , TT dA 



+ 



2! c/v?2 



+ 



3 rfM 



</3=0 



3! dip^ 



+ ... 



(/3=0 



</3=0 



2A^2 ^y,2 



+ 



d^A 



</3=0 



6Ar3 ^y,3 



+ ... (20) 



</3=0 



Comparing (!20|) with (flQll . if the lowest non- vanishing derivative at = is d^A/d(p^ for 
some integer j we have yh = j and a = {1 — j) / j . In order to find j, instead of using the 
comphcated solutions of (ITH]) . it is more elucidative to take consecutive derivatives of (ITSi) 
with respect to ip. From the first derivative of (ITS]) we deduce, with help of (fT^ - ffT7l) . at 

<p = : 



(Ai-A3)'(Ai-A4f(A3-A4)'Airf 







(21) 



(p=0 



Given that Ai = A2 at (^9 = 0, using = (rf^ — ^) in ([S]) we arrive at {d"^ — 1) (c^ — 1) = 
0. Since d and c are non negative numbers, see ([3]), and c = 1 corresponds to T — t- 00 
while d = 1 is the spin-1/2 Ising model for which a = —1/2 is known exactly , we assume 
henceforth 



Xfy^C^d'- 



d^ 



(22) 



Equations (ETj) and (122]) tell us that 



dA 

dip 



, if Ai = A2 7^ Aj , z = 3, 4 and A3 7^ A4 



(23) 



ip=0 



With help of (IT1D-([I7D and (ES]), the second derivative of ([H]) at = furnishes: 



From 



(Ai-A3)'(Ai-A4)'(A3-A4)'Aj2ci 



1 



Af 

c 



d^A 



d^^ 



1/3=0 



+A?[(Ai-A3)(Ai-A4)]} =0 



and (EH) we have : 



(24) 



<PA 

dip"^ 



7^ , if Ai = A2 7^ Aj , i = 3, 4 and A3 7^ A4 



i/)=0 



(25) 
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In summary, from fl23p and fl25p we conclude that yh = 2 a.s long as we have neither a triple 
degeneracy of the transfer matrix eigenvalues nor two double degeneracies. Therefore, for 
the one-dimensional ANNNI model we show on general grounds that a = —1/2, see also 
|17j . except for the two mentioned special cases where a different critical behavior may 
appear in principle. 

In the latter cases we have no information from (12T]) and f l24p . For triple degeneracy 
(Ai = A2 = A3) the coefficients in ([5]) must satisfy: 

18 — 14 010203 + 3 0103 — 0303 + 4 = 200(802 — 303] (26) 
9 [12 {aj - 01O2O3) + 3 oio^ - O2O3] + 32 o^ = (27) 

Working out f l26|) and f l27|) we arrive at : 

d^(8rf^-9)(l-c^)2 + 108c^(l-d^)2 = , (28) 



, 2 8^4-7 , , 

A = 29 

8(i4-9 ^ ^ 

Since both c and d are nonnegative real numbers, it is clear from (!28|) that we can only have 
triple degeneracy if < c/ < (9/8)1/^. The condition ([28]) is a second degree polynomial 
on c^. Thus, there are only two possibilities for the temperature as a function of d: 



c±{d) 



,3/2 , . 



{1-ld^y' ±{d''-i] 
(1 T(rf^-i) 



(30) 



Notice, in agreement with fl28|) . that c„((i) = l/c+{d). In figure 1 we plot both c±{d'^). 
They coalesce into c_ = c+ = 1 (T — )■ 00) at = 1. The function c+{d) diverges at 
d = (9/8)^''^ while c^{d) vanishes at that point. By inserting c^{d) in the exact solutions 
of the cubic equation flTSl) and expanding the results about = we obtain: 



2d^^ -40d^ + 117d^-81 9 (rf^-1) (2^4-9) (2d4_3) 



jd'-l) {3645 + 2d^ [-1944 + rf^ (297 + 64d^)]} 
+ ^Y6d^^ ^ (31) 
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2 (8 c^^ - 7) 32i{d^-l) V2#-3 
8d^-9 ^ 9 (8 - 9) 





(32) 



2 (8 rf^ - 7) _ 32 i (rf^ - 1) V2d^-3 
8#-9 9 (8 # - 9) 



8 (rf^ - 1) , 
27 



+ 




(33) 



The solutions A2 and ^3 are interchanged under the symmetry (f — )■ —(f of ( fT8l) while 
is invariant. 

It turns out that only A2{f) and A^^ip) satisfy, at (p = 0, the triple degeneracy 
condition (129|) . Indeed, substituting c = c+{d) and A = Ai{ip) in ([5]) one obtains the four 
transfer matrix eigenvalues Aq, = Xai'p) , a = 1,2,3,4 for each solution A = Ai{{p) , i = 
1, 2, 3. We have checked numerically for several values of d in the range < d < (9/8)^/^ 
that A = Ai{ip) only leads to double degeneracy Ai = A2 at ^9 = 0. Besides, it is such 
that |Ai((/9)| = |A2(v5)| < |Aj(v2)| , z = 3,4 in the neighborhood of ip = 0. So, Ai{0) does 
not correspond to a true edge singularity. On the other hand, the function A^lip), though 
it leads to triple degeneracy Ai = A2 = A3 at = 0, it is such that |Ai(v9)| = |A2(v5)| is 
not the largest absolute value in the vicinity of <y9 = 0. Thus, we do not have partition 
function zeros for A = A3 (</)). For A = ^2(^5) we have checked that Ai = A2 = A3 at 
ip = and more importantly |Ai(v9)| = |A2(v5)| > |Aj(v9)| , i = 3,4 which confirms that we 
do have Yang-Lee zeros approaching the YLES A2{0) for A = A2{'p). 

From the above discussion and (l32il we conclude that for the fine-tuning c = c^{d) 
we assure a different critical behavior with j = 3 = i/h {cr = —2/3) for the density of 
Yang-Lee zeros at the YLES. 

Regarding the second possibility c = c_((i) = l/c+(c?) it is possible to show that in 
this case we have the usual result j = 2 = yh {cr = —1/2). Indeed, it can be checked 
analytically that the cubic equation ( |T8l) is symmetric under c — )■ 1/c. However, there is 
no such symmetry in ([5]). If we plug c = c_{d) and Ai{(p) in ([5]), it turns out that Ai = A2 
at ip = and |Ai(v9)| = |A2(<y5)| > |Ai(v5)| , i = 3,4 in the neighborhood of = 0. So we 
have a true YLES with, see (EH), j = 2 = yh {a = —1/2). The other functions A2{p>) and 
Aslip) lead to Ai = A2 = A3 at = but |Ai(v9)| = |A2(v2)| is not the largest absolute 
value about ip = 0. So, we do not have Yang-Lee zeros in those cases. Analogously, the 
case of two double degeneracies Ai = A2 7^ A3 = A4 leads only to a = —1/2. 
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We see from f l29p . which gives the location of the YLES in the triply degenerated case, 
that we can only have a = —2/3 either for zeros lying on the unit circle (—2 < A < 2), 
with < d < 1 (AFM coupling < 0), or on the negative real axis {A < —2) which 
requires 1 < d < (9/8)^/"^ (FM coupling K > 0). In the first case < c+(d) < 1 (FM 
coupling J > 0) while in the second one c > 1 (AFM coupling J < 0). Therefore, the 
unusual critical behavior a = —2/3 only occurs for couplings J and K of opposite nature. 

As a final remark we note from (132|) that j ^ 3 if d^ = 3/2, which coincides precisely, 
using c = c+{d), with the quadruple degeneracy of the eigenvalues Ai = A2 = A3 = A4. 
However, in this case c+ becomes complex and it will be neglected here. Anyway, this is 
an indication that different values for a are associated with multiple degeneracies of the 
transfer matrix eigenvalues. 

In the next section our analytic results are confirmed by numerical calculations of the 
Yang-Lee zeros. 



4 Numerical results 

Comparing the FSS relation f|T9|) for two rings of sizes and Na+i we derive a numerical 
estimate for i/h 



yh = - 



In 



N, 



a+l 



In 



(34) 



Where either the imaginary or the real part of A'Ui(A^) can be used. In the case of 
triple degeneracy, the YLES uf{oo) = are known exactly by inverting the relation 
ue + 1/ue = Ae where Ae is given in fl29|) for each value of d. Since Ae E ^ it follows 
that 



[ui 



So we choose ue = without loss of generality. 



We also consider another finite size scaling relation 
zeros close to the YLES: 



l30] for the linear density of 



p{L) = C2 L^^'^ = C2Ny'^~^ 

where C2 is a constant independent on the number of the spins while A^ = 
Analogous to we can derive from (!35|l : 



(35) 
= L. 



yh = i + 



In- 



A^, 



a+l 



In 



(36) 
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The scahng exponents obtained from and will be called respectively and yf^. 
More specifically, we use y^'^'^ and y^'^"^ according to the use of real or imaginary parts 
of Aui{N). We stress that fl34l) and fl36|) furnish independent numerical estimates for 
yh, since p{N) = 1/ (A^|ui — U2I) depends upon the first and second closest zeros, ui , U2 
to the YLES while yj^ depends only upon the first one. Later on, we will extrapolate 
the finite size results fl51|) and fl5B]) for — )• 00 via BST (Burlish-Stoer) extrapolation 
algorithm [HI [32]. 

The partition function zeros for a ring with N sites (spins) are obtained numerically 
with help of the software Mathematica from an analytic expression for Z^. Even for 
one-dimensional spin models there are no analytic expressions for the Yang-Lee zeros in 
general. In order to save computer time, instead of using the analytic solution for given 
in in terms of the transfer matrix eigenvalues or in terms of the trace of powers of the 
transfer matrix as in Tr T^, we use an alternativ^ exact expression derived in pO] for any 
spin model which can be solved via a finite transfer matrix. Namely, since Aj , z = 1, 2, 3, 4 
are solutions of the secular equation P4(A) = A"^ — a^X^ + a2A^ — aiA + oq = 0, we have 
shown, formula (11) of [20j, that can be identified with 

Z^ = -N {In [g^P^il/g)] = -N [in (l - a,g + a^g^ - a,g^ + aog')]^^ , (37) 

where g is an arbitrary real variable (power counting parameter) and [f{g)]gN stands for 
the coefficient of the term of power g^ in the Taylor series of f{g) about g = For the 
lowest powers = 1,2,3 the reader can easily check, with help of (fT^ -(fT71) at = 0, 
that fl37|) indeed reproduces the transfer matrix solution 

At each value of d, the expression c+((i), see fl30l) . furnishes the corresponding fine- 
tuned temperature for triple degeneracy,. We can also invert c+{d) and obtain d for 
each given value of c. If, for instance, we choose c = 0.5, the inversion of c+{d) leads to 
d ~ 0.8680. We have displayed in figures 2-4 the YLZ and half of the corresponding YLES 
(with positive imaginary part) for d 0.8680 and c = 0.48 , 0.50 , 0.52. It turns out that 
the triple degeneracy point is a turning point after which each edge bifurcates into two 
new ones. Right above the triple degeneracy point (c > 0.50) we have checked (not shown 
here) that at the endpoint of each of the two new edges, figure 4(b), the critical exponent 
is the usual one yh = 2 (a = —1/2) while right before (c < 0.50) we have a crossover 
behavior fiowing from y^ = 2 to y^ = ?> a,s we approach c = 0.50 from below. 

■^The alternative formula (|37|) has a diagrammatic interpretation as a connected Feynman diagram of 
a zero-dimensional Gaussian field theory |20] . 
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At the triple degeneracy point we have made a detailed analysis of the scaling behavior 
of the zeros in the neighborhood of the YLES which is located, see fl29p at d ^ 0.8680, 
at Ae ~ 1.1031 {ue ~ 0.5516 ± 0.8341 i). In this case all Yang-Lee zeros lie on the unit 
circle, see figure 3. 

The log-log fits in figures 5(a) and 5(b) confirm the FSS relations fll9p and fl35p . They 
furnish the estimates yj^ = 2.9960 (using the real pat of the zeros) and = 2.9866. 

In table ([1]) we present the sequences yh{Na) obtained from formulae (IMI) and (l36l) . 
In the last line we have extrapolated our finite size results — )■ oo by using the BST 
[311[32] algorithm with a; = 1. This algorithm approximates the original sequence yh{Na) 
by another sequence of ratios of polynomials with faster convergence. The BST approach 
depends upon the real free parameter oo: yh{N) = yh{oc) + ^ + -^h + ••• where Ai, A2, ■ ■ ■ 
are iV-independent constants. We plot the extrapolated quantity yh{oo) for 0.1 < w < 3.0 
altogether with their error bars in figure 6. The error bar corresponds to twice the differ- 
ence between the values of yh obtained at the step before the last one in the extrapolating 
sequence. In table ([T]) we have chosen u = 1 because it provides a more stable result, 
i.e., ^ = 0. Clearly from table ([1]) and figures 5 and 6 we have a result very close to 
yh = S {a = —2/3) at the triple degeneracy point. We have also checked numerically 
other couples of values for (c, d) satisfying the triple degeneracy condition c = c+{d). The 
BST extrapolated results for yh are very similar as well as their error bars. Some caution 
is needed when the edges are nearly horizontal (vertical) lines. In those cases the smallest 
error bars for yj^ are obtained by the use of the real (imaginary) part of the first zero 
respectively. 



5 Conclusion 

Usually, in one- dimensional spin models, the linear density of partition function zeros 
diverges with a critical exponent a = — 1/2 at edge singularities. The universality of a 
is known for a long time [6l [TJ [12] and checked explicitly in D = 1 in several models, 
see e.g. [131 El ESI EHl EI]- However, in the works [HI [191 120] one has found another 
critical behavior (cr = —2/3). The models investigated in [TSl [T^l [20] have three-state 
per site and only nearest-neighbor interaction. Here we have shown that a = —2/3 also 
appears in the one-dimensional spin-1/2 ANNNI model which contains a next-to-nearest- 
neighbor interaction and only two states per site. Our results support the universality of 
cr = —2/3. As in [HI [20], the triple degeneracy of the transfer matrix (TM) eigenvalues is 
necessary to evade the well known result a = —1/2. Such condition requires a fine-tuning 
of the couplings of the model which explains why the authors of [17] have only found 
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cr = —1/2 for the same model treated here. So, rather than the number of states per site, 
the important point is the dimension of the TM and the number of the free parameters 
of the model to be fine-tuned. 

The above argument signalizes that the same phenomenon might occur in higher- 
dimensional spin models under special circumstances, since for D > 1 the number of 
eigenvalues of the TM increases with the size of the lattice. In particular, one could 
speculate that this phenomenon might be behind the sudden drop from a = —0.15(2) 
down to cr = —0.365 as reported in [11], similarly to the drop from a = — 1/2 to cr = —2/3. 
In [H] one obtains the linear density of Yang-Lee zeros for the two-dimensional Ising 
model, above the critical temperature, indirectly from a function that fits the experimental 
magnetization data from a sample of FeCl2. The 2D-Ising model works as a prototype 
for FeCl2 in some temperature range. As shown already in [6] the discontinuity of the 
magnetization across the curve of zeros furnishes their density. So one has indirect access 
to a experimentally. 

Another interesting point is that the triple degeneracy condition c = c^{d) as given in 
( 130|) defines a transition point between two different loci of Yang-Lee zeros. For c < c^{d) 
we have an arc of the unit circle with two edges while for c > c^{d) each edge bifurcates 
into two new edges with some fraction of zeros leaving the unit circle, see figure 2 and 
figure 4. Our figure 4 is similar to figure 2-c of p/Tj. We only have a = —2/3 at c = c+{d). 

At last, we remark that the subtle breakdown of the permutation symmetry between 
the two largest eigenvalues (Ai ^ A2) is the key point in finding the unusual critical 
behavior with yh = 3{a = —2/3). 
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c 




0.2 0.4 0.6 0.8 1 

Figure 1: The solid (dashed) hne represents c+ (c_) as given in formula fl30|) 



Im[u] 




(a) (b) 

Figure 2: (Color online) Yang-Lee zeros for = 60 spins at d ^ 0.8680 and c = 0.48. In 
figure (a) we have all zeros while figure (b) displays only the zeros on the first quadrant 
with the corresponding Yang-Lee edge singularity (green (light gray) dot). The solid line 
stands for ImI = 1 
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Im[u] 




(a) (b) 

Figure 3: (Color online) Yang-Lee zeros and half of the Yang-Lee edge singularities (green 
(light gray) dot) for = 60 spins at d ^ 0.8680 and c = 0.5 (triple degeneracy point) 
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ln[Re(Au)] 




(a) Straight line: InRe(Au) = 2.022 - (b) Straight line: In p = -5.341 + 1.987 In 
2.996 In TV with Xe = 9-30 x 10"^ with Xp = 9.03 x 10"^ 

Figure 5: Log-log fits of ([19]) and ([35]) at (i ^ 0.8680 and c = 0.5 with 160 < < 240 
spins 







E,Re 
Vh 


E.Im 

Vh 


160 


2.984155688412 


2.995280492398 


2.995273778849 


170 


2.985057758429 


2.995549136960 


2.995543509484 


180 


2.985862192536 


2.995788825553 


2.995784061898 


190 


2.986584090647 


2.996004003019 


2.995999935011 


200 


2.987235580790 


2.996198248276 


2.996194746764 


210 


2.987826518204 


2.996374475494 


2.996371439968 


220 


2.988364995418 


2.996535081652 


2.996532432944 


230 


2.988857720521 


2.996682056444 


2.996679731514 


oo 


3.00000000000(1) 


3.00000000000(1) 


3.00000000000(3) 



Table 1: Finite size results for the Yang-Lee zeros, where yf^ are obtained from (1361) 
while y^'^^ and y^'^"^ come from (IMj) using the real and imaginary parts of the zeros 
respectively. The data of this table have been obtained from rings with 160 < N < 240 
spins at c = 0.5 and d ^ 0.8680. The last row is the N ^ oo extrapolation via BST 
algorithm with u = 1. 
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3.0002 

3 

2.9998 
^ 2.9996 
2.9994 
2.9992 
2 . 999 



♦♦♦♦ ♦ ♦ 



3.001 



3 * ♦ , 



ci^ 2.999 
2.998 
2.997 



0.5 1 1.5 2 2.5 3 



1.5 1 1.5 2 2.5 3 



(a) 



(b) 



ure 6: BST extrapolation for 0.1 < u < 3.0 of yj^ using the real part of the zeros 
ure 6(a) ) and (figure 6(b) ) at ci ~ 0.8680 and c = 0.5. 
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